%matplotlib notebook
from ipywidgets import interact, SelectionSlider, IntSlider, FloatSlider
import matplotlib.pyplot as plt
from matplotlib import animation
import numpy as np
import scipy.signal
from IPython.display import YouTubeVideo, HTML, Audio
from bokeh.layouts import column, row
from bokeh.models import CustomJS, ColumnDataSource, Slider
from bokeh.plotting import Figure, show, output_notebook
output_notebook()
t = np.arange(-5, 5, step=0.0001)
s = lambda t, sigma=0.5 : np.exp(-0.5*t**2/sigma**2)
Fs = Slider(start=0.2, end=5, value=1, step=.01, title="Frecuencia de muestreo")
t_dirac = np.arange(-5, 5, step=1/Fs.value)
create_figure = lambda title : Figure(plot_width=280, plot_height=280, title=title,
toolbar_location="below", x_range=(-5, 5))
p1 = create_figure('Señal Gaussiana')
p2 = create_figure('Peineta de Dirac')
p3 = create_figure('Señal muestreada')
p1.line(t, s(t), line_width=3)
for p in [p1, p2, p3]:
p.xaxis[0].axis_label = 'Tiempo [s]'
source = ColumnDataSource(data=dict(t=t_dirac,
tops=np.ones_like(t_dirac),
sd=s(t_dirac)
))
p2.vbar(x='t', top='tops', source=source, width=0.05)
p3.scatter('t', 'sd', source=source)
callback = CustomJS(args=dict(source=source, Fs=Fs), code="""
var t_dirac = [];
var tops = [];
var sd = []
for (let i = -5.0; i <= 5.0; i+=1/Fs.value) {
t_dirac.push(i);
tops.push(1);
sd.push(Math.exp(-0.5*Math.pow(i, 2)/Math.pow(0.5, 2)));
}
source.data['t'] = t_dirac;
source.data['tops'] = tops;
source.data['sd'] = sd;
source.change.emit();
""")
Fs.js_on_change('value', callback)
show(column(Fs, row(p1, p2, p3)))
%%capture
fig, ax = plt.subplots(2, figsize=(7, 4), sharex=True, tight_layout=True)
t = np.arange(-4, 4, step=1e-4)
def my_signal(t, a=0, T=0.5):
s = np.zeros_like(t)
s[np.absolute(t-a)<T] = 1
return s
# 2 segundos de espacio entre los dientes de la peinte
def peineta(t):
s = np.zeros_like(t)
s[::20000] = 1
return s
conv_s = np.convolve(my_signal(t), peineta(t), mode='same')
def update(a = 0):
ax[0].cla(); ax[1].cla()
p1, p2 = my_signal(t, 0.1*a - 4), peineta(t)
ax[0].plot(t, p1, label='señal');
ax[0].plot(t, p2, label='peineta');
ax[0].legend()
ax[1].plot(t, conv_s);
ax[1].set_xlabel('Tiempo [s]');
ax[1].scatter(0.1*a -4, np.sum(p1*p2), s=100, c='k')
return ()
anim = animation.FuncAnimation(fig, update, frames=80, interval=100, blit=True)
HTML(anim.to_html5_video())
x = np.arange(-5, 5, step=0.001)
sigma = Slider(start=0.1, end=2, value=1/np.sqrt(2*np.pi), step=.01, title="sigma")
source = ColumnDataSource(data=dict(x=x,
gt=np.exp(-0.5*x**2/sigma.value**2),
gf=np.exp(-2*(np.pi*x*sigma.value)**2)*np.sqrt(2*np.pi)*sigma.value
))
create_figure = lambda title : Figure(plot_width=350, plot_height=280, title=title,
toolbar_location="below", x_range=(-5, 5))
p1 = create_figure('Dominio de tiempo', )
p2 = create_figure('Dominio de frecuencia')
p1.line('x', 'gt', source=source, line_width=3)
p2.line('x', 'gf', source=source, line_width=3)
callback = CustomJS(args=dict(source=source, sigma=sigma), code="""
var x = source.data['x'];
var gt = source.data['gt'];
var gf = source.data['gf']
for (var i = 0; i < x.length; i++) {
gt[i] = Math.exp(-0.5*Math.pow(x[i]/sigma.value, 2));
gf[i] = Math.sqrt(Math.PI*2)*Math.exp(-2*Math.pow(Math.PI*x[i]*sigma.value, 2))*sigma.value;
}
source.change.emit();
""")
sigma.js_on_change('value', callback)
show(column(sigma, row(p1, p2)))
def gaussiana(f, sigma):
return np.sqrt(2*np.pi*sigma**2)*np.exp(-2*(np.pi*f*sigma)**2)
def espectro_discreto(f, sigma):
S = np.zeros_like(f)
for m in range(-20, 20):
S += gaussiana(f - Fs*m, sigma)
return S
def ventana(f, Fs):
SW = np.zeros_like(f)
SW[np.absolute(f) < Fs/2] = 1
return SW
Fs = 2.
sigma = 1.
f = np.arange(-3*Fs, 3*Fs, step=1e-3)
p1 = Figure(plot_width=600, plot_height=280, toolbar_location="below")
p1.line(f, gaussiana(f, sigma), color='green', alpha=0.75,
line_width=3, legend_label='Espectro original')
p1.line(f, espectro_discreto(f, sigma), line_width=3, alpha=0.75,
legend_label='Espectro discreto')
p1.line(f, ventana(f, Fs), line_width=3, alpha=0.75, color='black',
line_dash='dashed', legend_label='Ventana cuadrada')
p1.xaxis[0].axis_label = 'Frecuencia [Hz]'
show(p1)
Fs = 2.
sigma = 0.4
f = np.arange(-3*Fs, 3*Fs, step=1e-3)
p1 = Figure(plot_width=600, plot_height=280, toolbar_location="below")
p1.line(f, gaussiana(f, sigma), color='green', alpha=0.75,
line_width=3, legend_label='Espectro original')
p1.line(f, espectro_discreto(f, sigma), line_width=3, alpha=0.75,
legend_label='Espectro discreto')
p1.line(f, ventana(f, Fs), line_width=3, alpha=0.75, color='black',
line_dash='dashed', legend_label='Ventana cuadrada')
p1.xaxis[0].axis_label = 'Frecuencia [Hz]'
show(p1)
Fs = 0.75
sigma = 1.
f = np.arange(-3*Fs, 3*Fs, step=1e-3)
p1 = Figure(plot_width=600, plot_height=280, toolbar_location="below")
p1.line(f, gaussiana(f, sigma), color='green', alpha=0.75,
line_width=3, legend_label='Espectro original')
p1.line(f, espectro_discreto(f, sigma), line_width=3, alpha=0.75,
legend_label='Espectro discreto')
p1.line(f, ventana(f, Fs), line_width=3, alpha=0.75, color='black',
line_dash='dashed', legend_label='Ventana cuadrada')
p1.xaxis[0].axis_label = 'Frecuencia [Hz]'
show(p1)
import scipy.fft as sfft
f0 = 1.23 # Frecuencia de la sinusoide
T = 100 # Largo temporal
Fs = 2 # Frecuencia de muestreo
ts = np.arange(0, T, step=1/Fs)
signal = lambda t, f : np.cos(2.0*np.pi*f*t)
SA = sfft.fftshift(np.absolute(sfft.fft(signal(ts, f0))))
freq = sfft.fftshift(sfft.fftfreq(n=len(ts), d=1/Fs))
p1 = Figure(plot_width=600, plot_height=280, toolbar_location="below")
p1.line(freq, SA, line_width=3)
p1.xaxis[0].axis_label = 'Frecuencia [Hz]'
show(p1)
f = np.arange(-3*Fs, 3*Fs, step=1e-3)
def espectro_coseno(f, f0):
S = np.zeros_like(f)
S[np.isclose(f, -f0)] = 1
S[np.isclose(f, f0)] = 1
return S
def espectro_discreto(f, f0):
S = np.zeros_like(f)
for m in range(-20, 20):
S += espectro_coseno(f - Fs*m, f0)
return S
def ventana(f, Fs):
SW = np.zeros_like(f)
SW[np.absolute(f) < Fs/2] = 1
return SW
p1 = Figure(plot_width=600, plot_height=280, toolbar_location="below")
p1.line(f, espectro_coseno(f, f0), color='green', alpha=0.75,
line_width=4, legend_label='Espectro original')
p1.line(f, espectro_discreto(f, f0), line_width=3, alpha=0.75,
legend_label='Espectro discreto')
p1.line(f, ventana(f, Fs), line_width=3, alpha=0.75, color='black',
line_dash='dashed', legend_label='Ventana cuadrada')
p1.xaxis[0].axis_label = 'Frecuencia [Hz]'
show(p1)
T = 5
t = np.arange(0, T, step=1e-4)
ts = np.arange(0, T, step=1/Fs)
signal = lambda t, f : np.cos(2.0*np.pi*f*t)
p1 = Figure(plot_width=600, plot_height=280, toolbar_location="below")
p1.line(t, signal(t, f0), color='green', alpha=0.75,
line_width=3, legend_label='Señal original')
p1.scatter(ts, signal(ts, f0), size=10, legend_label='Señal muestreada')
p1.line(t, signal(t, Fs - f0), line_width=3, alpha=0.75, color='black',
line_dash='dashed', legend_label='Alias de la señal')
p1.xaxis[0].axis_label = 'Tiempo [s]'
show(p1)